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Abstract 

Charged particle production in neutral current deep inelastic scattering (DIS) has 
been studied using the ZEUS detector. The evolution of the mean multiplicities, 
scaled momenta and transverse momenta in Q 2 and x for 10 < Q 2 < 5120 GeV 2 
and x > 6 x 10~ 4 has been investigated in the current and target fragmentation 
regions of the Breit frame. Distributions in the target region, using HERA data 
for the first time, are compared to distributions in the current region. Predictions 
based on MLLA and LPHD are inconsistent with the data. 
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1 Introduction 



This paper reports the results of a study of the properties of the hadronic final state 
in neutral current positron-proton deep inelastic scattering (DIS). The fragmentation of 
the struck quark in DIS is compared to that of the quarks produced in e + e~ annihila- 
tion, allowing the concept of universality of fragmentation in different processes to be 
tested. A comparison is also made of the fragmentation associated with the outgoing 
struck quark to that in the target region; the latter is associated with the initial state 
radiation from the incoming proton. The data are compared to analytical Quantum 
Chromodynamic (QCD) predictions for the momenta spectra for both the current and 
target region. 

The event kinematics of DIS are determined by the negative square of the four-momentum 
of the virtual exchanged boson, Q 2 = —q 2 , and the Bjorken scaling variable, x = 
Q 2 /2P-q, where P is the four-momentum of the proton. In the Quark Parton Model 
(QPM), the interacting quark from the proton carries four- momentum xP. The variable 
y, the fractional energy transfer to the proton in its rest frame, is related to x and 
Q 2 by y ~ Q 2 /xs, where y/s is the positron-proton centre of mass energy. The invari- 
ant mass, W, of the hadronic system is related to x, Q 2 and the proton mass, m p , by 
W 2 = Q 2 (l — x)/x + m 2 At fixed Q 2 , the W 2 behaviour reflects the x dependence. 
A natural frame in which to study the dynamics of the hadronic final state in DIS is the 
Breit frame [JjJ. In this frame, the exchanged virtual boson (7*) is completely space-like 
and has a four-momentum q = (0, 0,0, — Q = —2xP Bre%t ) = (E, p x , p y , p z ), where 
pBreit j g momentum of the proton in the Breit frame. The particles produced in the 
interaction can be assigned to one of two regions: the current region if their ^-momentum 
in the Breit frame is negative, and the target region if their z-momentum is positive (see 
Fig. [I]). The main advantage of this frame is that it gives a maximal separation of the 
incoming and outgoing partons in the QPM. In this model the maximum momentum a 
particle can have in the current region is Q/2, while in the target region the maximum 
is Q{\ — x)/2x. In the Breit frame, unlike the hadronic centre of mass (7*2?) frame, the 
two regions are asymmetric, particularly at low x, where the target region occupies most 
of the available phase space. 

The measurements presented here extend the previous studies of fragmentation in the 
Breit frame performed at HERA |2]44|]. Increased statistics coming from an integrated 
luminosity of 38 pb _1 lead to a significant improvement in the precision of the scaled 
momentum distributions of the charged particles in the current fragmentation region. 
A subset of data, corresponding to 6.4 pb _1 , has also been used to measure a small 
part of the target fragmentation region (the dark shaded region indicated in Fig. [l]); 
the study of the target region is limited by acceptance and systematic effects rather 
than by statistics. The scaled and the transverse momentum distributions of charged 
particles in the hadronic final state are measured in the current and target regions 
of the Breit frame as a function of x and Q 2 in the ranges 10 < Q 2 < 5120 GeV 2 
and x > 6 x 10~ 4 . Comparisons are made with Monte Carlo models, QCD analytical 
calculations and e + e~ data. 
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2 QCD Models 



The current region in the ep Breit frame is analogous to a single hemisphere of e + e~ 
annihilation. In e + e~ — > qq annihilation the two quarks are produced with equal and 
opposite momenta, ±^/s ee /2. The fragmentation of these quarks can be compared to that 
of the quark struck from the proton; this quark has an outgoing momentum —Q/2 in the 
Breit frame. In the direction of this struck quark the scaled momentum spectra of the 
particles, expressed in terms of x p = 2p Breit /Q, are expected to have a dependence 
on Q similar to that observed in e + e~ annihilation at energy ^Js ee = Q, with no x 
dependence. The effects of higher order processes not present in e + e~ annihilation are 
discussed in Ref. ||]. 

Results from e + e _ annihilation support the need for coherence effects in perturbative 
QCD [S— 13 1 . The phenomenon of coherence is a natural consequence of the quantum 
mechanical treatment of the parton cascade. Long wavelength gluons are unable to 
resolve individual colour charges of partons within the parton cascade. This has the 
effect that the available phase space for soft gluon emissions is reduced to an angular- 
ordered region, due to destructive interference. This leads to a number of important 
differences in the properties of the partonic final state relative to the incoherent case. 
The most notable of these are a slower rise in the multiplicity of partons with increasing 
initial parton energy, and the modification of the parton logarithmic momentum spectra 
to a skewed Gaussian form, often referred to as the "hump-backed" plateau [12]. The 
parton level predictions in practice depend on two free parameters, a running strong 
coupling, governed by a QCD scale A, and an energy cut-off, Qq, below which the 
parton evolution is truncated. In this case A is an effective scale parameter and is not 
to be identified with the standard QCD scale, e.g. A^g. In particular, predictions can 
be made at Qq = A yielding the so-called limiting spectrum. 

Within the framework of the modified leading log approximation (MLLA) there are 
predictions of how the higher order moments of the parton momentum spectra should 
evolve with the energy scale [14, 15 [. The MLLA calculations are made at the parton 
level. The hypothesis of local parton hadron duality (LPHD) fll6|| , which relates the 
observed hadron distributions to the calculated parton distributions via a constant of 
proportionality, is used in conjunction with the parton predictions of the MLLA to 
allow the calculation to be directly compared to data. There is uncertainty about the 
energy scale at which the applicability of LPHD breaks down, in which case the parton 
predictions cannot be compared to the hadron distributions. 

DIS at low x allows a study of fragmentation in the target region following the initial 
scattering off a sea quark (or antiquark). The description based on MLLA [ IT | is shown 
schematically in Fig. ^, where the quark box at the top of the gluon ladder represents the 
scattered sea quark plus its antiquark partner. There are various contributions to these 
calculations. Contribution C, the top leg of the quark box, corresponds to fragmentation 
in the current region. Three further contributions (Tl, T2 and T3), which are sources 
of soft gluons, are considered in these analytical calculations to be associated with the 
target region. It is predicted [|l7j that the contribution Tl behaves in the same way as 
the current quark C and so should have no x dependence. The contribution T2 is due 
to the colour field between the remnant and the struck quark, and the contribution T3 
corresponds to the fragmentation of the rungs in the gluon ladder. Both T2 and T3 are 
predicted to have x and Q 2 dependences which differ from Tl. Both the Tl and T2 
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contributions have been calculated and give particles of momenta < Q/2. The collinear 
gluons T3, on the other hand, generally fragment to particles with momentum For 
values of the scaled momentum x p < 1.0, the region of phase space is analogous to the 
current region and has contributions mainly from Tl and T2. The parton momentum 
spectra predicted by MLLA, over a range of Q 2 and x, are shown in more detail in 
Ref. Q. In the target region these spectra are approximately Gaussian for x p < 1; they 
peak at a value of x p ~ 0.1 — 0.2 in the range of x and Q 2 measured in this paper, falling 
to a plateau region for 1 < x p < (1 — x)/x (the maximum value of x p in the target 
region). The acceptance of the ZEUS detector allows the study of the contributions 
from C, from Tl and from T2. 

Scaling violations are predicted in the fragmentation functions, which represent the 
probability for a parton to fragment into a particular hadron carrying a given fraction 
of the parton's energy. Fragmentation functions incorporate the long-distance non- 
perturbative physics of the hadronisation process in which the observed hadrons are 
formed from final state partons of the hard scattering process. Like parton densities, 
they cannot be calculated in perturbative QCD but can be evolved with the hard- 
process scale, using the DGLAP evolution |l8| equations, from a starting distribution at 
a defined energy scale; this starting distribution can be derived from a fit to data. If the 
fragmentation functions are combined with the cross sections for the inclusive production 
of each parton type in the given physical process, predictions can be made for scaling 
violations, expressed as the Q 2 evolution of the x p spectra of final state hadrons [19]. 
These scaling violations allow a measurement of the strong coupling constant, a s , and 
such studies have been performed at LEP |p0| , [2l]l by incorporating lower energy PETRA 
data. The NLO calculations (CYCLOPS) p!| of the scaled momentum distribution exist 
for DIS but as yet the appropriate fragmentation functions at different values of Aj^g 
do not exist to allow the extraction of a s from DIS data. 



3 Experimental Setup 

The data presented here were taken at the positron-proton collider HERA using the 
ZEUS detector. The 38 pb _1 of data corresponds to data taken in 1994-1996 and part 
of the 1997 data sample. The 1995 data alone, corresponding to 6.4 pb" 1 , was used to in- 
vestigate the target region. During the period 1994-1997 HERA operated with positrons 
of energy E e = 27.5 GeV and protons with energy 820 GeV. The ZEUS detector is a 
multipurpose detector. Of particular importance in this analysis are the central tracking 
detector (CTD) and the uranium-scintillator calorimeter (CAL). A detailed description 
of the ZEUS detector can be found in p3|,p?4||. 

Throughout this paper we use the standard ZEUS right-handed coordinate system, in 
which X = Y = Z = 0is the nominal interaction point, the positive Z-axis points in 
the direction of the proton beam (referred to as the forward direction) and the X-axis 
is horizontal, pointing towards the centre of HERA. 

The tracking system consists of a central tracking chamber (CTD) p5| in a 1.43 T 
solenoidal magnetic field. The CTD, which surrounds the beampipe, is a drift chamber 
consisting of 72 cylindrical layers, arranged in 9 superlayers. Superlayers with wires 
parallel to the beam axis alternate with those inclined at a small angle to give a stereo 
view. The single hit efficiency of the CTD is greater than 95%. The resolution of 
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the transverse momentum, p t (lab), in the laboratory frame for full-length tracks can 
be parametrised as a(p t (lab))/p t (lab) = 0.0058p t (lab) © 0.0065 © 0.0014/p t (lab), with 
Pt(l&b) in GeV. (This form is a more precise description of the transverse momentum 
resolution, particularly at low to medium pt(lab), than that given hitherto). 
Surrounding the solenoid is the uranium-scintillator calorimeter (CAL) which is 
divided into three parts: forward, barrel and rear covering the polar regions 2.6° to 
36.7°, 36.7° to 129.1° and 129.1° to 176.2°, respectively. The CAL covers 99.7% of the 
solid angle, with holes in the centres of the forward and rear calorimeters to accom- 
modate the HERA beampipe. Each of the calorimeter parts is subdivided into towers 
which are segmented longitudinally into electromagnetic (EMC) and hadronic (HAC) 
sections. These sections are further subdivided into cells each of which is read out by 
two photomultipliers. Prom test beam data, energy resolutions of <te/E = 0.18/a/ E for 
electrons and cte/E = 0.35/VE~ for hadrons (E in GeV) have been obtained. 
The ZEUS detector is almost hermetic, allowing the kinematic variables x and Q 2 to 
be reconstructed in a variety of ways using combinations of energies and angles of the 
positron and hadronic system. Variables calculated only from the measurements of 
the energy, E' e , and angle, 6> e , of the scattered positron are denoted with the subscript 
e, whilst those calculated from the hadronic system measurements, with the Jacquet 



Blondel method [27], are denoted by the subscript JB. Variables calculated by these 
approaches are used only in the event selection. In the double angle method p8[ , denoted 
by DA, the kinematic variables are determined using 6 e and the angle 7# (which is the 
direction of the struck quark in the QPM), defined from the hadronic final state. 
The DA method was used throughout this analysis for the calculation of the boosts 
and the kinematic variables because it is less sensitive to systematic uncertainties in the 
energy measurement than other methods. 

The triggering and online event selections were identical to those used for the measure- 



ment of the structure function F% [29|. The reconstructed tracks used in the charged par- 
ticle analyses are associated with the primary event vertex and have 7^ (lab) > 150 MeV 
and |r/(lab)| < 1.75, where r/(lab) is the pseudorapidity given by — ln(tan(0/2)) with 
6 being the polar angle of the measured track with respect to the proton direction in 
the lab. frame. This is the region of CTD acceptance where the detector response and 
systematics are best understood. 

Further selection criteria were applied both to ensure accurate reconstruction of the 
kinematic variables and to increase the purity of the sample by eliminating background 
from photoproduction processes: 

• E' e > 10 GeV, to achieve a high purity sample of DIS events; 

• Q 2 da — 10 GeV 2 , to further enhance the purity of the DIS sample; 

• Ve < 0.95, to reduce the photoproduction background; 

• Ujb > 0.04, to give sufficient accuracy for DA reconstruction of Q 2 and x; 

• 35 < 5 = (E — pz) < 60 GeV where the summation is over all calorimeter cells, 
to remove photoproduction events and events with large radiative corrections; 

• \X\ > 16 cm or \Y\ > 16 cm, where X and Y are the impact position of the 
positron on the CAL, to avoid the region directly adjacent to the rear beampipe; 



4 



• — 40 < Vertex < 50 cm, to reduce background events from non-ep collisions. 

The (x, Q 2 ) bins are listed in Table [l|. The sizes of the bins were chosen to give good 
statistics in each bin and to limit the migrations between bins [0]. There is negligible 
background from non-DIS events. 



4 Event Simulation 

Monte Carlo event simulation is used to correct for acceptance and resolution effects. 



The detector simulation is based on the GEANT 3.13 [3C] program and incorporates 
our best knowledge of the apparatus. 

To calculate the correction factors, neutral current DIS events were generated, via the 
DJANGO 6.24 program pi, using HERACLES 4.5.2 p2| which incorporates first order 



electroweak corrections. The QCD cascade was modelled with the colour dipole model, 
including the boson-gluon fusion process, using the ARIADNE 4.08 program. In this 
model coherence effects are implicitly included in the formalism of the parton cascade. 
The colour dipole model treats gluons emitted from quark-antiquark (diquark) pairs as 
radiation from a colour dipole between two partons. This results in partons that are not 
ordered in their transverse momenta. The program uses the Lund string fragmentation 
model [34 1 for the hadronisation phase, as implemented in JETSET 7.4 (J35jj . For the 
analysis of the 1995 data, two Monte Carlo samples were generated, 4.2 pb _1 with Q 2 > 6 
GeV 2 and 15.8 pb" 1 with Q 2 > 40 GeV 2 , using the GRV94 |36| parameterisation of the 
parton distribution functions. For the 1996 and 1997 data, a sample with Q 2 > 70 
GeV 2 , (with MRSA parton densities J37||), was generated, corresponding to 17.1 pb _1 , 
and a sample with Q 2 > 800 GeV 2 (with GRV94 parton densities), corresponding to 
53 pb" 1 . 

For the studies of the systematics for the 1995 data, two additional samples of events 
were generated (2.1 pb _1 with Q 2 > 6 GeV 2 and 9.1 pb^ 1 with Q 2 > 70 GeV 2 ) using 



the HERWIG 5.8d Monte Carlo program 38], where no electroweak radiative correc- 



tions were applied. In HERWIG, coherence effects in the QCD cascades are included by 
angular ordering of successive parton emissions and a clustering model is used for the 



hadronisation |ll],[3§. For the 1996 and 1997 data, HERWIG samples with Q 2 > 70 
GeV 2 were generated, corresponding to 9.0 pb -1 , and Q 2 > 800 GeV 2 (both with MRSA 
parton densities), corresponding to 60 pb -1 . Both the GRV94 and MRSA parametrisa- 
tions agree well with the HERA measurements of the proton structure function F2 in 
the (x, Q 2 ) range of this analysis ji(], 41 1. 



Another approach to modelling the parton cascade is included in the LEPTO 6.5.1 [42] 
program, which incorporates the LO a s matrix element matched to DGLAP parton 
showers (MEPS). This recent version of LEPTO incorporates the soft colour interaction 



(SCI) model [43] to describe HERA rapidity gap events. SCI produces changes to the 
usual string topologies in non-gap events causing the string to overlap itself and this 
results in an increase both of particle number and energy per unit of rapidity. 
The Linked Dipole Chain model, LDC 1.0, |44|] has also been investigated. In this model 
the parton shower evolution is based on a reformulation jig of the CCFM approach [46] 



which approximates the BFKL Wm prediction at low x and the DGLAP prediction in 



the high x limit. The parton density parametrisation of "set A" [44] was used, which 



fits data from HI and ZEUS. The DGLAP equation predicts strong ordering of the 
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parton transverse momenta while the BFKL equation relaxes this ordering but imposes 
strong ordering of the longitudinal momenta. Both the LEPTO and LDC programs 
use the Lund string fragmentation model. They were used to compare generator level 
calculations with our data. 

5 Correction Procedure 

The Monte Carlo event samples were used to determine the mean charged particle 
acceptance in the current region as a function of (x,Q 2 ). The chosen analysis intervals 
in (x, Q 2 ) correspond to regions of high acceptance (between 74 and 96%) in the current 
region of the Breit frame. The acceptance for the limited area of the target region under 
study (x p < 1.0) is lower; the ln(l/x p ) distributions have a good acceptance around 
their peak positions (70-90%) but it falls to about 50% at lower values of ln(l/x p ) for 
the (x, Q 2 ) bins 1-4 defined in Table |l|. Due to the low acceptance and large systematic 
uncertainties for ln(l/x p ) < 1.0 this region is not studied in (x,Q 2 ) bins 5-8 and no 
studies in the target region are made beyond bin 8. There is good acceptance (>90%) 
for both regions under study for pt > 1.0 GeV, where pt is the transverse momentum with 
respect to the virtual photon direction in the Breit frame. However the acceptance falls 
below 50% for pt < 0.5 GeV. These values are well understood in terms of geometrical 
acceptances. 



Bin no. 


x range 


Q 2 (GeV 2 ) range 


1 


(6.0- 12.0)10~ 4 


10 - 20 


2 


(1.2 - 2.4)10~ a 


10 - 20 


3 




20-40 


4 




40-80 


5 


(2.4- 10.0)10" 3 


20-40 


6 




40 - 80 


7 




80 - 160 


8 




160 - 320 


9 


(1.0 - 5.0)10" 2 


160 - 320 


10 




320 - 640 


11 




640 - 1280 


12 


0.025 -0.15 


1280 - 2560 


13 


0.05 - 0.25 


2560 - 5120 



Table 1: The (x,Q 2 ) analysis bins. 



About 7% of the tracks generated in the current region migrate to the target region. 
Migrations into the current region from the target fragmentation region are typically 
less than 5% of the tracks assigned to the current region for Q 2 > 320 GeV 2 . For 10 < 
Q 2 < 320 GeV 2 these migrations are on average 12%, reaching 25% for Q 2 < 40 GeV 2 
and low values of y where the hadronic activity is low and the measurement of "fn is 
subject to systematic problems leading to a worse x resolution and hence an uncertainty 
in the boost vector to the Breit frame. 

The correction procedure is based on the detailed Monte Carlo simulation of the ZEUS 
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detector with the event generators described in the previous section. Since the ARI- 
ADNE model gives the best overall description of our observed energy flow |1 it is 
used for the standard corrections to the distributions. 

The data are corrected for trigger and event selection cuts; event migration between 
(x,Q 2 ) intervals; QED radiative effects; track reconstruction efficiency; track selection 
cuts in pt(lab) and r/(lab); track migration between the current and target regions; and 
for the products of A and K s decays which are assigned to the primary vertex. 
Correction factors were obtained from the Monte Carlo simulation by comparing the 
generated distributions, excluding decay products of A and Kg, with the reconstructed 
distributions after the detector and trigger simulations. The same reconstruction, selec- 
tion and analysis were used for the Monte Carlo simulated events as for the data. The 
correction factors, F{x p ), were calculated for each x p bin using a bin-by-bin correction: 

F(x p ) ' 





-ZVgen \ dx pJ gen / N obs \^P/ ohs 

where N gen (N Q ^ S ) is the number of generated (observed) Monte Carlo events in each 
(x, Q 2 ) interval and n is the number of charged particles (tracks) in the current or target 
region in the corresponding x p and (x, Q 2 ) interval. A similar correction procedure was 
applied for the other variables. The bin sizes of the distributions were chosen to be 
commensurate with the measurement resolution. In the current region, the overall 
correction factors are greater than unity and typically < 1.3. In the target region these 
correction factors are larger but, in the region that we measure, they are typically < 1.5 
for bins 1-4 and < 2.0 for bins 5-8. 



6 Systematic Checks 

The systematic uncertainties in the measurement can be divided into three types: uncer- 
tainties due to event reconstruction and selection, to track selection, and to the Monte 
Carlo model used. The systematic checks were as follows: 

• The cut on y e < 0.95 was changed to y e < 0.8. 

• The cut on 35 < S(E - p 2 ) < 60 GeV was changed to 40 < S(E - pj < 60 GeV. 

• The tracking cuts on |r/(lab)| < 1.75 and pt(lab) > 150 MeV were tightened to 
|ry(lab)| < 1.5 and pt(lab) > 200 MeV; the cuts were also removed. 

• Instead of requiring vertex fitted tracks, all reconstructed tracks that passed 
through superlayers one and three of the CTD were accepted. 

• The data were corrected using a different hadronisation model, namely HERWIG, 
in place of ARIADNE. 

With the exception of the change in the model from ARIADNE to HERWIG, and the 
use of non-vertex tracks, all the systematic effects were small, i.e. within two standard 
deviations of the statistical errors. Of the two major systematic uncertainties, the 
hadronisation model change was dominant. 
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Current fragmentation region 

The use of non- vertex tracks resulted in an overall increase in the single particle densities 
of 5% to 15% and was fairly flat across the x p range. The use of HERWIG to unfold the 
data gave rise to systematic shifts as large as 15%. For Q 2 < 80 GeV 2 the tendency of 
the correction was to lower the single particle density values at low x p and to increase the 
values at higher x p . For pt in the range < p% < 0.5 GeV, the systematic uncertainties 
were about 10%. They reduced with increasing p t to about 5%. 



Target fragmentation region 

Systematic effects due to the different hadronisation models were largest (as high as 
50%) in the ln(l/x p ) distributions at high x p . For x p < 0.3 they were typically 10%, 
increasing to 30% at larger x p . The tracking systematic from non- vertex tracks was 
largest at low x p in the lowest Q 2 bin where it is 7%, but otherwise was of the order of 
2%. The pt distribution was little influenced by the model used, as would be expected 
from the good acceptance, at all but the very lowest transverse momenta. The mean pt 
showed a model dependence in the target region of at most 20% which was due to the 
reduced acceptance for low p t tracks at high x p . In general, unfolding with HERWIG 
resulted in a higher value of the normalised single particle densities. 



7 Results 

7.1 Current fragmentation region 

Figure |3| shows the ln(l/x p ) distributions for charged particles in the current fragmen- 
tation region of the Breit frame for different bins of (x,Q 2 ). These distributions are 
approximately Gaussian in shape with the mean charged multiplicity given by the in- 
tegral of the distributions. As Q 2 increases, the multiplicity increases and, in addition, 
the peak of the distribution moves to larger values of m(l/x p ). The moments of the 
ln(l/x p ) distributions have been investigated up to the 4th order; the mean (/), width 
(w), skewness (s) and kurtosis (k) were extracted from each distribution by fitting a 
distorted Gaussian of the following form: 

1 ' ,a oce^lh-U-^ + ^ + ^ + ^-ksA, (1) 



a tot dln(l/x p ) *V8 2 4 V ' 6 24 

where 5 = (ln(l/a; p ) — l)/w, over a range of ±1.5 units (for Q 2 < 160 GeV 2 ) or ±2 units 
(for Q 2 > 160 GeV 2 ) in ln(l/x p ) around the mean. The equation is motivated by the 



expression used for the MLLA predictions of the spectra 14]. The smooth curves in 
Fig. ^ result from the fit of equation (|lj) to the data; they represent the data well. 
Figure || shows the moments of the ln(l/x p ) spectra as a function of Q 2 . It is evident 
that the mean and width increase with increasing Q 2 , whereas the skewness and kurtosis 
decrease. Similar fits performed on e + e~ data show a reasonable agreement with 
our results, consistent with the universality of fragmentation for this distribution. 
The data are compared to the MLLA predictions of Ref. [15|, using a value of A = 



175 MeV, for different values of Qq. A comparison is also made with the predictions of 
Ref. [[l4j] for the limiting spectrum (Qq = A). The MLLA predictions of the limiting 
spectrum in Ref. jLq] describe the mean well. However both of the MLLA calculations 
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predict a negative skewness which tends towards zero with increasing Q 2 in the case of 
the limiting spectra. This is contrary to the measurements. The qualitative description 
of the behaviour of the skewness with Q 2 can be achieved for a truncated cascade 
(Qo > A), but a consistent description of the mean, width, skewness and kurtosis 
cannot be achieved. A range of A values was investigated and no single value of A 
gave a consistent description of all the moments. 

We conclude that the MLLA predictions, assuming LPHD, do not describe the data. 
We note however that a moments analysis has been performed [j50|, taking into account 
the limitations of the massless assumptions of the MLLA predictions, and yields good 
agreement between the limiting case of the MLLA jn| and e + e~ data over a large energy 
range, 3.0 < y / s^ < 133.0 GeV. A discussion of phase space effects on the ln(l/x p ) 



distributions is given in ref. [51|. These phase space effects can resemble MLLA. 
In Fig. H the evolution of the moments with Q 2 (same DIS data as Fig. are compared 
with the predictions of various Monte Carlo models. Both ARIADNE and LEPTO 
(with SCI) give a reasonable description of the data, while HERWIG fails to predict the 
observed Q 2 variation. This is particularly noticeable for skewness and kurtosis. The 
discontinuities in the HERWIG prediction arise from a strong x dependence in bins of 
overlapping Q 2 . Such an x dependence is not observed in the data. It may be noted that, 
due to the choice of the maximum scale of the parton shower evolution, there are fewer 
gluons radiated in HERWIG than in the other generators; this could possibly account 
for the poor agreement of HERWIG with our measurement. All Monte Carlo programs 
have been compared using the default values of their parameters. The LEPTO model 
without SCI (not shown) describes the data better than does the default version. 
The inclusive charged particle distribution, 1/crtot da/dx p , in the current fragmentation 
region of the Breit frame is shown in bins of x p and Q 2 in Fig. ||. The fall-off as Q 2 
increases for x p > 0.3 (corresponding to the production of more particles with a smaller 
fractional momentum) is indicative of scaling violations in the fragmentation function. 
The distributions rise with Q 2 for x p < 0.1 and are discussed in more detail below. The 
data are compared to e + e~ data []52| (divided by two to account for the production of a 
qq pair) at Q 2 = s ee . For the higher Q 2 values shown there is a good agreement between 
the measurements in the current region of the Breit frame in DIS and the e + e~ results; 
this again supports the universality of fragmentation. The fall-off observed in the ZEUS 
data at low x p and low Q 2 is greater than that observed in e + e~ data at SPEAR ]53fl ; 
this can be attributed to processes not present in e + e~ (e.g. scattering off a sea quark 
and/or boson gluon fusion (BGF)) which depopulate the current region 



A kinematic correction has recently been suggested jgjj to the NLO calculation []22j] of 
the inclusive charged particle distribution which has the form: 

1 



(2) 



where m e ff is an effective mass to account for the massless assumption used in the 
fragmentation functions. It is expected to lie in the range 0.1 GeV < m e g- < 1.0 GeV. 
The x p data are compared to the CYCLOPS NLO QCD calculation incorporating this 
correction in Fig. [?]. This calculation convolutes the fragmentation function of each 
type of parton with the cross sections for their production. It combines a full next-to- 
leading order matrix element with the MRSA' parton densities (with Aqcd = 230 MeV) 
and NLO fragmentation functions derived from fits to e + e~ data pq]. The kinematic 
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correction allows a more legitimate theoretical comparison to lower Q 2 and x p than 
was possible in our earlier publication ||. The bands represent the uncertainty in the 
predictions by taking the extreme cases of m e fr = 0.1 GeV and m e fj = 1.0 GeV. These 
uncertainties are large at low Q 2 and low x p , becoming smaller as Q 2 and x p increase. 
Within these theoretical uncertainties there is good agreement throughout the selected 
kinematic range. The kinematic correction describes the general trend of the data but 
it is not possible to achieve a good x 2 fit for m e fl over the whole x p and Q 2 range. The 
uncertainties introduced by these additional processes restrict to high Q 2 and high x p 
the kinematic range that may be used to extract a s from the observed scaling violations. 
The pt distributions, l/<Ttot da/dp 2 , are shown in Fig. |8| for x p < 1.0. The distributions 
show an exponential fall off at low pt although it is evident that a high-p^ tail develops 
with increasing Q 2 . These high pt tails contribute at most 15% of the cross section. 
The pt distributions, from the 1995 data, in the first 8 bins of (x,Q 2 ) are shown in 

Fig. H as closed data points. The straight lines are exponential fits, exp(—b^Jp 2 + m 2 ), 
to the low pt interval 0.2 - 1.0 GeV, where is the mass of the pion. They yield 
slopes of b ~ 5 — 6 GeV" 1 . The values of b show little Q 2 dependence. For bin 8 the line 
extending out to higher pt (pt > 1.2 GeV) is a fit to the empirical power law formula [57] 



A X (1 +Pt/pto)~ m - There are strong correlations between A,p t o and m. Consequently 
Pto has been fixed at 0.75 GeV, a value consistent with the fit with all variables free and 
that used by HI in the analysis of their photoproduction data [55]. With this pto, the 
parameterisation fits the data well (x 2 /NdF = 5.1/19) and gives m = 5.8 ± 0.4 ± 0.1. 
These tails are slightly higher if the x p < 1.0 cut is removed (open points in the figure). 
Particles with x p > 1.0 occur due to hard QCD processes, such as BGF and QCD 
Compton. 



bin no. 


<Q 2 > 
(GeV 2 ) 


< n > ± stat ± syst 


1 


14.0 


1.13 ±0.01 ±0.05 


2 


14.1 


1.18 ±0.01 ±0.04 


3 


27.9 


1.70 ±0.01 ± 0.07 


4 


55.3 


2.27 ±0.01 ±0.07 


5 


28.0 


1.81 ±0.01 ±0.06 


6 


55.9 


2.44 ±0.01 ±0.14 


7 


110. 


3.00 ±0.01 ± 0.23 


8 


216. 


3.77 ±0.02 ± 0.26 


9 


221. 


3.98 ±0.02 ± 0.37 


10 


443. 


4.59 ±0.03 ± 0.40 


11 


863. 


5.26 ±0.05 ± 0.39 


12 


1766. 


6.01 ±0.05 ± 0.46 


13 


3507. 


7.10 ±0.11 ±0.69 



Table 2: Mean charged multiplicity in the current fragmentation region. 



Figure [H] shows the mean charged multiplicity in the current fragmentation region. The 
results, for each of the (x,Q 2 ) bins, are listed in Table [2]. The multiplicity increases by 
about a factor of six over the measured Q 2 range. Also shown in Fig. 1C are results from 
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e + e~ annihilation experiments |5^| (scaled down by a factor of 2) and results from fixed 
target DIS data |30| at similar Q 2 (Q 2 < 30 GeV 2 ) to the ZEUS data but corresponding 
to an x range about two orders of magnitude higher. For Q 2 ,>80 GeV 2 there is reasonable 
agreement between the results from e + e~ data and ZEUS, again consistent with the 
universality of fragmentation. At lower Q 2 the multiplicities measured by ZEUS are 
lower than those found in the e + e~ data and the fixed target DIS data. Similar results 
have recently been observed by the NOMAD collaboration |6l]]. In this Q 2 region there 
is a negligible contribution from charmed quarks so that the difference must originate 
from the depopulation of the current region due to the prevalence of boson-gluon fusion 



processes in this low (x, Q ) region ||,|T3],|62|. 



Figure 11 displays the same ZEUS data as in Fig. [H] but now compared to various Monte 
Carlo models. Both ARIADNE and HERWIG, with default settings, describe well the 
variation of the multiplicity with Q 2 . LEPTO with SCI, while describing the data at 
low Q 2 , simulates the Q 2 evolution incorrectly which leads to an overestimation of the 
multiplicity at high Q 2 . This overestimation of the data by LEPTO can be partially 
rectified, as can be seen in Fig. |ll|(dash-dotted line), by removing SCI from the model. 

7.2 Target fragmentation region 

The distributions of charged particles in x p and transverse momentum, pt, in the target 
region of the Breit frame are studied as a function of x and Q 2 , in the range 6 x 10 -4 < 
x < 1 x 1CT 2 and 10 < Q 2 < 320 GeV 2 . This analysis, which uses the 1995 data, is 
restricted to values of the scaled momenta x p < 1.0 so that similar phase space regions 
for the target and current can be compared. Thus the accepted part of the target region 
corresponds to the contributions Tl and T2 in Fig. |2] and is only a small part of the 
complete phase space, as depicted in Fig. [j]. This is also the kinematic region of the 
target fragmentation that has a reasonable acceptance. The corrected data distributions 
with combined statistical and systematic errors are shown in Figs. 12-i~6|. 
The distributions in ln(l/x p ) are shown for both the target and current regions in 



Fig. 12. The fitted curves shown are two-piece normal distributions |5^] to guide the 
eye. In contrast to the current region, the target region distribution does not fall to 
zero as ln(l/x p ) tends to zero. Although the magnitude of the single particle density 
at the peak position of the current region distribution grows by a factor of about three 
over the Q 2 range shown, the single particle density of the target distribution, at the x p 
value corresponding to the peak of the current distribution (contribution C is equivalent 
to contribution Tl in Fig. ^), depends less strongly on Q 2 and increases by only about 
30%. In addition the ln(l/x p ) distribution shows no significant dependence on x when 
Q 2 is kept constant. In the target region the peak position of the ln(l/x p ) distribution 
increases more rapidly with Q 2 than in the current region; this is consistent with the 
behaviour expected from cylindrical phase space. The approximate Gaussian distribu- 
tion of the MLLA predictions peaking at ln(l/x p ) ~ 1.5 — 2.5 || is not observed. We 
conclude that the target distributions are inconsistent with the MLLA predictions when 
used in conjunction with LPHD. 

The ln(l/x p ) distributions in the target fragmentation region are compared to Monte 



Carlo models in Fig. [TJ. The ARIADNE and Linked Dipole Chain (LDC) models 
describe the data well in the measured (x, Q 2 ) bins. The two Monte Carlo models 
based on DGLAP parton evolution techniques, LEPTO and HERWIG, fail to describe 
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bin no. 


< Q 2 > (GeV 2 ) 


current < n > 


target < n > 


1 


14.0 


1.13 ±0.01±0.05 


4.95 ±0.01±0.2 


2 


14.1 


1.18 ±0.01±0.04 


4.94 ±0.01±0.4 


3 


27.9 


1.70 ±0.01±0.07 


6.11 ±0.02±0.6 


4 


55.3 


2.27 ±0.01±0.07 


7.36 ±0.03±0.3 



Table 3: The mean charged multiplicities in the current and target regions for the 
(x,Q 2 ) analysis bins, in the range < x p < 1. The first error is statistical and the 
second is systematic. 



the data. The LEPTO Monte Carlo with SCI describes the data at low Q 2 but the 
dependence on Q 2 within the model is incorrect and discrepancies are observed at large 
Q 2 . The HERWIG Monte Carlo gives a poor description of the data in all (x, Q 2 ) bins. 
The LEPTO generator without SCI (not shown) gives a good description of the data. 
The mean multiplicity in the target region in the range < x p < 1, shown in Table ||, is 
larger than in the current region, by about a factor of four in the lowest (x, Q 2 ) bins. The 
target region multiplicity increases with Q 2 , but more slowly than that in the current 
region presumably due to the additional fragmentation terms in the target region shown 
in Fig. ^. Only the first four (x,Q 2 ) bins are studied as they have reasonable acceptance 
over the whole of < x p < 1. Fig. [l4] shows the ratio of the charged multiplicities in 
the target and the current regions as a function of Q 2 . The ratio falls as Q 2 increases. 
Also shown are comparisons with Monte Carlo models; ARIADNE, LDC and LEPTO 
without SCI all describe the trend of the data. HERWIG, though reproducing the Q 2 
dependence, fails to predict the magnitude whilst LEPTO with SCI fails to describe the 
Q 2 dependence of the ratio. 

The pt distributions, 1 / atotdcr / dp 2 , are shown in Fig. 15 and the same fits have been 
performed on these distributions as on the current region pt distributions of Fig. 0. 



The fit of the exponential, exp(—byp 2 + m 2 ), gives slopes of b ~ 5 — 6 GeV -1 . In 
a similar manner to the distributions in the current region, the values of b exhibit 
little Q 2 dependence, but the distributions develop a high-p t tail with increasing Q 2 . 
The line plotted for the Q 2 interval 160-320 GeV 2 is a fit to the power law formula 
Ax(l+p t /p t0 )~ m for p t greater than 1.2 GeV. This fits the data well (x 2 /NdF = 1.9/19) 
with pto fixed at 0.75 GeV and m = 5.7 ± 0.6 ± 0.1. These high pt tails contain at most 
15% of the cross section. The values of m in the current and target regions agree within 
errors and are smaller than that found for the ZEUS j5?J and HI |58| photoproduction 
data. This is consistent with what would be expected from the point-like nature of the 
exchanged photon in DIS. 

To compare the general characteristics of the transverse momentum distributions in the 
target and current regions, the mean pt versus x p is shown in Fig. 16. The mean pt at 
large x p is higher in the current region than in the target region and shows a stronger 
Q 2 dependence than the target region. In the current region the mean pt rises with x p 
and reaches a maximum as x p tends to —1. In contrast, in the target region the mean 
Pt tends to a constant value of about 0.6 GeV. Thus the target region p t distribution 
has the characteristics of pt-limited phase space with only a small dependence on Q 2 . 
The mean pt vs x p distribution in the target region shows no significant dependence on 
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x when Q 2 is kept constant. 

Also shown in Fig. are the Monte Carlo predictions of the ARIADNE and HERWIG 
models. The ARIADNE generator gives a good description of the data and is very 
similar to the predictions of LEPTO and LDC Monte Carlo models. The HERWIG 
generator gives a less satisfactory description of the data. The discrepancy in the target 
region can partially be explained by the lack of intrinsic transverse momentum of the 
incoming struck quark in the default parameters of HERWIG. 

The apparent contradiction between the similarity of Figs. |9| and 15 and the difference 
between the mean pt in the current and target regions in Fig. 1£ may be understood 
in terms of the correlation between x p and pt- In the current region this correlation 
is strong with the high x p ,pt region corresponding to low multiplicity. As the high pt 
particles have a strong Q 2 dependence, this is reflected in the Q 2 dependence of the 
mean pt at high x p . In contrast, in the target region the correlation between x p and 
Pt is small. This results in a lower mean pt at large x p and a mean pt substantially 
independent of x p and Q 2 as x p tends to 1. 



8 Summary 

Charged particle distributions have been studied in the Breit frame in DIS over a wide 
range of Q 2 . The distributions in scaled momentum, x p , and transverse momentum, 
Pt, have been measured for the first time in the target region of the Breit frame for 
1.2 x 1CT 3 < x < 1 x 1CT 2 and 10 < Q 2 < 320 GeV 2 . For scaled momenta in the interval 
< x p < 1 the mean target region charged track multiplicity is found to be larger than 
that measured in the current region; there is no significant x-dependence at fixed Q 2 . 
The transverse momentum distributions for both the current and target fragmentation 
regions exhibit similar properties. A tail at large pt develops as Q 2 increases. The mean 
transverse momentum as a function of x p has a weaker dependence on Q 2 in the target 
region than the current region. Whereas in the current region the mean pt increases 
approximately linearly with x p , the mean pt in the target region tends to a constant 
value with increasing x p , consistent with cylindrical phase space. 

The HERWIG model is unable to describe the Q 2 dependence of the ln(l/x p ) distri- 
butions in the target fragmentation region. In contrast the colour dipole model as 
implemented in the ARIADNE program, LEPTO generator without SCI and the LDC 
Monte Carlo, based on a reformulation of the CCFM evolution, all adequately describe 
the data. 

In the current region, the results show clear evidence for scaling violations in scaled 
momenta as a function of Q 2 and support the hypothesis of the coherent nature of 
QCD cascades. The data are well described by NLO calculations. The comparison 
of our results in the current region of the Breit frame with e + e~ data at Q 2 = s ee 
for Q 2 > 80 GeV 2 shows good agreement. The moments of the ln(l/x p ) spectra in 
the current region exhibit the same energy scale behaviour as those observed in e + e~ 
data. The observed charged particle spectra are consistent with the universality of quark 
fragmentation in e + e" and DIS at high Q 2 . The moments cannot be described by the 
MLLA calculations together with LPHD. 

The target region ln(l/x p ) distribution shows a weaker Q 2 dependence than the corre- 
sponding current region distribution. In particular, the magnitude of the single particle 
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density of the target distribution, at the x p value corresponding to the peak of the cur- 
rent distribution, increases by about 30%, in contrast to a threefold increase for the 
current region in the Q 2 -range considered here. The MLLA predictions for the target 
region, in conjunction with LPHD, fail to describe the data. 
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CURRENT TARGET 

Figure 1: The phase space diagram for DIS in the Breit frame, pi denotes the lon- 
gitudinal momentum axis, referred to as the z— direction, and pt denotes transverse 
momentum. The dark shaded region indicates the part of the target region under study 
in this paper. 
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Figure 2: A schematic of DIS scattering at low x within the MLLA framework. Quark 
C represents the struck sea quark in the current fragmentation region. Tl is the other 
half of the quark box which is in the target region. T2 is the t-channel gluon exchange 
and T3 the rungs of the gluon ladder. 
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Figure 3: The charged particle distributions of \ / a tot da / dhi(l / x p ) in the current frag- 
mentation region as a function of ln(l/x p ) for different (x,Q 2 ) bins. Only statistical 
errors are shown. The full line is the skew Gaussian fit; the arrows indicate the fit 
range. 
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Figure 4: Evolution of the mean, width, skewness and kurtosis of the ln(l/x p ) distri- 
bution in the current fragmentation region with Q 2 . Data from e + e _ and ep are shown 
together with the MLLA predictions of Dokshitzer et al fl5| (the full line is Qo = A, 
the dashed Qq = 2A, and the dotted Qq = 3A) and the limiting spectrum predictions 
of Fong and Webber Jl4| (dash-dotted line where available.) The overlapping points are 
different x ranges in the same Q 2 range. The inner error bars are the statistical error 
and the outer error bars are the systematic and statistical errors added in quadrature. 
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Figure 5: Evolution of the mean, width, skewness and kurtosis of the ln(l/x p ) dis- 
tribution in the current fragmentation region with Q 2 . The overlapping points are for 
different x ranges in the same Q 2 range. The inner error bars are the statistical error 
and the outer error bars are the systematic and statistical errors added in quadrature. 
The lines are the predictions from the Monte Carlo generators ARIADNE (full), LEPTO 
with SCI (dashed) and HERWIG (dotted). The shaded region represents the statistical 
error from the fits to the Monte Carlo simulations. The LEPTO model without SCI 
resembles the predictions of ARIADNE. 
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Figure 6: The inclusive charged particle distribution, 1/atot da/dx p , in the current 
fragmentation region of the Breit frame. The inner error bar is the statistical and the 
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Figure 8: The transverse momentum distributions in the current fragmentation region 
for the 1994-1997 data (x p < 1.0) for different regions of x and Q 2 . The outer error bars 
are the statistical errors; the inner error bars are the sum of statistical and systematic 
errors added in quadrature. 
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Figure 9: The transverse momentum distributions in the current fragmentation region 
for the 1995 data for different regions of Q 2 and x. The inner error bars are the statistical 
errors; the outer error bars are the sum of statistical and systematic errors added in 
quadrature. The lines are the fits discussed in the text. The closed data points are for 
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Figure 10: The mean charged multiplicity for the current region. The error bars are the 
sum of the statistical and systematic errors in quadrature. The open points represent 
data from e + e~ experiments divided by two to account for q and q production (also 
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Also shown are fixed target DIS data. 



26 



ZEUS 1994-97 



E 

<D 

i_ 
O 

o 

c 
o 




12 



10 - 



8 - 



6 - 



4 - 



2 - 



ARIADNE 

LEPTO (with SCI) 

HERWIG 

LEPTO (no SCI) 

• ZEUS 94-97 (ep) 




1 10 1 2 1 3 1 4 

Q 2 (GeV 2 ) 
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Figure 12: The corrected ln(l/x p ) distributions for the target and current regions for 
the 1995 data. Fitted two-piece normal distributions are shown to guide the eye. The 
heavy line corresponds to the target region, the light line to the current region. The 
error bars are the sum of the statistical and systematic errors in quadrature. 
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Figure 13: The corrected ln(l/x p ) distributions for the target fragmentation region for 
the 1995 data compared to Monte Carlo models: ARIADNE, LEPTO with SCI, LDC 
and HERWIG. The error bars are the sum of the statistical and systematic errors in 
quadrature. The LEPTO model without SCI resembles the predictions of ARIADNE. 
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are the sum of the statistical and systematic errors in quadrature. The discontinuities, 
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Figure 15: The transverse momentum distributions in the target fragmentation region 
for the 1995 data. The inner error bars are the statistical errors; the outer error bars 
are the sum of statistical and systematic errors added in quadrature. The lines are the 
fits discussed in the text. 
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Figure 16: The corrected mean transverse momentum versus scaled momentum dis- 
tributions for the 1995 data. The inner error bars are the statistical errors; the outer 
error bars are the sum of statistical and systematic errors added in quadrature. The 
full line indicates the ARIADNE Monte Carlo prediction and the dotted the HERWIG 
prediction. 



32 



